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00 The new methodology of adding QCD NLO corrections in the initial state Monte Carlo parton 

t ] shower (hard process part) is presented using process of the heavy boson production at the LHC 

as an example. Despite the simplified model of the process, presented numerical results prove 
that the basic concept of the new methodology works correctly in the numerical environment of 

,__( the Monte Carlo parton shower event generator The presented method is an alternative to the 

^ well established methods, MC@NLO and POWHEG. Refinements of the new method with better 

• ^ computer CPU time efficiency are also discussed. 
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1. Introduction 

The Large Hadron Collider (LHC) at CERN provides rich harvest of experimental data. The 
proper understanding and interpretation of these data, possibly leading to discovery of new phe- 
nomena, requires perfect mastering of the "trivial" effects due to the multiple emissions of soft and 
coUinear gluons and quarks. Perturbative Quantum Chromodynamics (pQCD) [1, 2, 3], supple- 
mented with clever modelling of the low energy nonperturbative effects, is an indispensable tool 
for disentangling the Standard Model physics component in the data. This work presents part of 
the global effort of improving quality of the pQCD calculations for LHC experiments. 

Most of the results presented here are described in refs. [4] and [5]. Although this work elab- 
orates on the improved method of the pQCD calculation combining NLO-corrected hard process 
and LO parton shower Monte Carlo (MC), it should be regarded as the first step towards NNLO- 
corrected hard process combined with the NLO parton shower MC [6] . 

2. Basic LO parton shower MC 

The multigluon distribution of the single initial state ladder, which is a building block of our 
parton shower MC, is represented by the integrand of the "exclusive/unintegrated PDF", which in 
the LO approximation is the following: 

D{t,x) = J dxo dZ 5^=.voZ dQ{io,XQ) G{t,to - lnxQ\Z), 

GitMZ) = e-'^ I f n / d^m %>4_. ^^fe)) (2.1) 

where evolution kernel hP{z) = 5 (1 +z^), evolution time isfo = ln(^o/A) and the "erkonal" phase 
space integration element is d^S{k) = p = T^^^ ^"^^ ~ ^k^. We use rapidities 
E,i = i In ^ I j^j^ in the hadron beam rest frame (Rh), and T],- = ^ In ^ Irfhp defined in hai^d process 

rest frame (RFHP). They are related by = In ^ — T],. Rapidity ordering is now t = (^^ax > > 
■■■ > > ^i-i > ■ ■ ■ > ^0 = tQ, where to = ^0 = \n{qo/mi,) — InxQ. The direction of the z axis in the 
RFHP is pointing out towards the hadron momentum. A lightcone variable of the emitted gluon is 
defined as a,- = ^ and of the emitter parton (quark) as Xj = xq — L}=o (after / emissions). We 
also use fractions Zi = Xj/xj^i. The Sudakov formfactor Sp comes from the "unitarity" condition' 
Jq dZ G{t,to\Z) = 1, which is also instrumental in the Markovian MC implementation used to 
obtain D{t,x) at any value oft> to. 

The initial distribution do{qo,xo) related to experiment, to previous steps in the MC ladder, or 
to PDF in the standard MS system is not essential for the following discussion, we only note that the 
unitarity condition provides baryon number conservation sum rule dxD{t,x) = Jq dxQ do{to,xo). 

For testing our new method of correcting hard process to the NLO level we use the follow- 
ing simplified MC parton shower, implementing the DY process with two ladders and the hard 

'The usual cutoff 1 — z < e regularizing the IR singularity is implicit. 
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Figure 1: The upper plot shows the LO distribution of rj^ ~ ^ \n{xf /xb) from the CMC LO parton shower 
(purple) and from the strictly collinear formula (green). The lower plot shows the ratio of the two. 
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ni+ii2 J — 

:dZ2{P- Y kj;qi,q2)^{sxFXB,e)W^^''. 



In the LO approximation W^^*^ = 1. Rapidity ^ is translated into 77 - the center of mass system 
rapidity, in the forward part (F) of the phase space as = In ^ — T],-, T]of > T]i > S, and in the 
backward (B) part as i^, = — In-^ + T],-, S > T], > tjob- The rapidity boundary between the two 
hemispheres S = is used, until a more sophisticated version related to rapidity of the produced 
VK/Z is introduced. 

Analytical integration of eq. (2.2) results in the standard factorization formula (W^^'^ = I) 



(7o = / dXf dXB DF{t,Xf) DB{t,XB) OB{sXfXB) 

Jo 



(2.3) 



The distributions Df{t,XF) = {do®Gf){t,Xf) and DB{t,XB) = {do®GB){t,XB) are obtained from 
separate Markovian LO Monte Carlo runs. The above LO formula is exact, and can be tested with 
an arbitrary numerical precision. 

Figure I represents a "calibration benchmark" for the overall normalization at the LO level. 
We show there the properly normalized distribution of the variable rj^ = \ \n{xF /xb), which in the 



2Following ref. [4], we adopt dr^iP^quqi) = 8^^\P-qx -qi)'^'^- 
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coUinear limit approximates the rapidity of W boson. The distribution in the upper plot of Fig. 1 , 
representing eq. (2.3), is obtained using the general purpose MC program FOAM [7]. The coUinear 
PDF D{t,x) there has been obtained from a separate high statistics run (10^° events) of a Markovian 
MC (MMC), creating D{t,x) in a form of the 2-dimensional look-up table^. The other distribution 
in the upper plot of Fig. 1 represents eq. (2.2) in LO approximation. It comes from the full scale 
MC generation (with four-momenta conservation). The MC run with 10^ events was used. The 
constrained MC (CMC) technique of ref. [10] is used here because of the narrow Breit-Wigner 
peak due to a heavy boson propagator^. Two CMC modules and FOAM are combined into one MC 
generating gluon emissions and the W boson production. FOAM is taking care of the generation of 
the variables Xf,XB,XfQ,XBo and the sharp Breit-Wigner peak ms = sxfXB, then two CMC modules 
are initialized and generate the gluon four-momenta . They are mapped into , following the 
prescription defined in ref. [4], such that the overall energy-momentum conservation is achieved. 
Figure 1 demonstrates a very good numerical agreement between da/drj^ from our full scale LO 
parton shower MC of eq. (2.2) and the simple formula of eq. (2.3), to within 0.5%, as seen from 
the ratio of the two results in the lower part of the figure. 



3. Introducing NLO corrections to hard process 

The NLO corrections to hard process are imposed on top of the LO distributions of eq. (2.2) 
using a single "monolithic" weight W^l^^ defined exactly as in ref. [4] : 

jeFP{ZFj) dOB{s,d)/dQ. jeBPyZBj) dOB{s,d)/dQ. 
the NLO soft-i- virtual correction is Ay+s = (i^^ ~ ?)' '•^^ ^^^^ correction reads: 



l + jl-a-Pf dOB l + il-a-PfdOB ^ 

-^«>'^ 2 da,^''^^-^'^<^ 2 do^}''^^- 



(3.2) 



The above is the exact ME of the quark-antiquark annihilation into a heavy vector boson with 
additional single real gluon emission^. The LO component, which is already included in the LO 
MC, is subtracted here. The variable s = sxfXB = {qi + ^2)^ is the effective mass squared of the 
heavy vector boson. The definition of angle 6 in the LO component is rather arbitrary. We define 
it in the rest frame of the heavy boson, where qi +q2 = 0, as an angle between the decay lepton 
momentum qi and the difference of momenta of the incoming quark and antiquark 6 = Z(qi, pof — 
Pob)- On the other hand the two angles in the NLO ME are defined quite unambiguously as df = 
Z{q[,—pQB) and Ob = Z{qi,pof). In the above we only need directions of the pof and pqb vectors, 
which aie the same as the directions of the hadron beams. The lightcone variables aj and j3y of the 

■'This MMC run solves the LO DGLAP equation using the MC method, as in refs. [8, 9]. 

backward evolution algorithm of ref. [11] could be also used here. 
^We employ here the compact representation of ref. [12], which has also been used in POWHEG [13], 
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Figure 2: The upper plot shows (— ) NLO correction (alone) to the distribution of Tj^ = ^\n{xf/xB) from 
the parton shower (purple) and from the the strictly collinear formula (green). Their ratio is plotted below. 

emitted gluon are defined in the F and B parts of the phase space as follows^: 

aj = \- ZFj, = ccj , for j e F, 

Pj = l-ZBj, a, = i8;^>-2(i.--), for JgB. 

Again, the exact phase space integration of eq. (2.2) including Wj^^'^ of eq. (3.1) is feasible, 
and the resulting compact expression for the total cross section is obtained [4] : 

Oi = J^^ dXf dXBdzDp{t,Xf) DB{t,XB) OB{sZXFXB){5-=l{l+As+v)+C2r{z)}, (3-3) 

where C2,-(z) = ^ [-^{l-z)]. 

3.1 Numerical test of NLO correction 

Figure 2 represents a principal proof of concept of our new methodology for implementing 
the NLO corrections to the hard process in the parton shower MC. The plotted NLO correction 
to the rj^ distribution^ comes from the parton shower MC with the NLO-coiTccted hard process 
according to eqs. (2.2) and (3.1). Additionally we also plot there result of a simple collinear formula 
of eq. (3.3), where two collinear PDFs are convoluted with the analytical coefficient function C2r(z) 
for the hard process. Both results coincide within the statistical enor, see their ratio in the lower 
part of Fig. 2. 

Technically, the inclusion of the NLO correction in our parton shower MC is rather straight- 
forward, and is obtained by including W^l^'^ weight of eq. (3.1). MC is providing both LO and 
NLO-corrected results in a single run with weighted events. The NLO weight is strongly peaked 
near W^^^ = 1, positive, and without long-range tails. Its distribution is shown in Fig. 3. 

*See ref. [4] for more explanations. 

^Extra minus sign introduced to facilitate visualization. 
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Figure 3: The distribution of the NLO weight W^^° of eq. (3.1) 



In all numerical results we have set ^y+s = 0> as it is completely unimportant for the presented 
analysis. The initial distributions Jo(^05-^o) defined in ref. [5]. 



4. Simplification of the method and comparison with other methodologies 




Figure 4: The inclusive distribution of gluons on the log Sudakov plane of rapidity t — ^ and v = ln(l — z) 
(LHS). Contributions from all gluons weighted with the component weight Wj^^" (RHS). 

Our new method for introducing NLO corrections in the hard process, proposed in ref. [4] and 
tested in ref. [5], is an alternative to the two well established MC@NLO [14] and POWHEG [15, 
16] methodologies. With MC numerical implementation at hand, let us elaborate on the differences 
with the above two techniques in particular with the POWHEG technique. We shall also see that it is 
possible to make our method more efficient in terms of CPU time consumption. This improvement 
is not so critical in the present case of NLO corrected hard process, but may be quite useful in the 
case of correcting evolution kernels to the NLO in the ladder parts of the MC [6]. 

The most important differences with the POWHEG and MC@NLO techniques are: 

• The summation over all emitted gluons, without deciding which gluon is the one involved in 
the NLO correction and which ones are merely "LO spectators" in the parton shower. 
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• The absence of ( 1 / ( 1 — z) ) + distributions in the real part of the NLO corrections (virtual+soft 
correction is kinematically independent). 

To explain more clearly how W^l^^ of eq. (3.1) is distributed over the multigluon phase space, we 
restrict now to single ladder (hemisphere) with a simplified weight: 



W, 



NLO 



MC 



NLO 



W 



NLO 



P{ZF})dOB{s,e)/d^' 



(4.1) 



In order to find out the phase space regions specific for NLO corrections we consider inclusive 
distributions of gluons on the Sudakov logarithmic plane of rapidity t, and variable v = ln(l — z). 
In the left hand side (LHS) of Fig. 4 we show gluons inclusive distribution in the LO approximation. 
The flat plateau there represents IR/collinear singularity^ ICp^dE, with the drop by factor 1/2 
towards z = 0, due to factor in the LO kernel. In the right hand side (RHS) of Fig. 4 we show 
contributions from all gluons weighted with the component weight^ —^ft/i^^o gq (-4 jj^g NLO 
contribution is concentrated in the area near the hard process rapidity t = <§max, which has to be 
true for the genuine NLO contribution^". The completeness of the phase space near this important 
region (z = 0, i^max) is critical for the completeness of the NLO corrections. Both POWHEG and 
MC@NLO use standard LO MCs which feature an empty "dead zone" in this phase space corner. 

Figure 4 suggests that the dominant contribution to W^^^ could be from the gluon with the 
maximum ln/:J ~ i^y -|-ln(l —Zj), which is closest to the hard process phase space corner. In the 
MC we may easily relabel generated gluons using new index K such that they are ordered in the 
variable Kjf = c^^f -|-ln(l — za'), 'Qf+i < with K = \ being the hardest one. 



LO gluon K=l 



LO gluon K>1 



0.1^- 

o.oa] 

0.064 





Figure 5: The inclusive LO distribution of Fig. 4 split into the hardest gluon (left) and the rest (right). 

Figure 5 demonstrates a split of the LO inclusive distribution of Fig. 4 into the ^ = 1 compo- 
nent and the rest K> \. The important point is that the ^ = 1 component reproduces the original 
complete distribution over the whole region where the NLO correction is non-negligible! This is 
exactly the observation on which POWHEG technique is built. According to the POWHEG au- 
thors, taking the K= \ component is sufficient to reproduce the complete NLO correction (modulo 
NNLO). 

^We use constant 

^We again insert a minus sign in order to facilitate visualization, 
'"it also vanishes towards the soft limit z ^ 1. 
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Figure 6: The original NLO correction from Y.j VVj^^*^ and its two hardest gluon components, from W^^'^ 



and W^^2^ as a function of x = YljZj- 



The above statement is checked numerically in Fig. 6, where we compare the NLO correction 
to the X = YljZj distribution from the complete sum Wj^^^ and from Wj^^f. As we see the K = I 
component saturates the complete sum very well, with the K = 2 component being negligible in 
the first approximation. 

We can therefore speed up the calculation by means of taking only the K =1 contribution. The 
price will be that the formula of eq. (3.3) will not be exact any more. Our method differs, however, 
from the POWHEG scheme, where the ^ = 1 gluon is generated separately in the first step, and 
other gluons are generated (by the LO parton shower MC) in the next step. That is easy for LO MC 
with /:^-ordering, while in case of the LO MC with angular-ordering POWHEG requires additional 
effort of generating the so called vetoed and truncated showers. In our method, there is no need for 
such vetoed/truncated showers in case of angular- ordering. 



ALL: p (t,v) = .°" ; Plateau at 2C %= 0.10 




Angular ordering: J=l 



Angular ordering: J=2 




0.08- 
0.06- 
0.04 
0.02 



Angular ordering: J=7 ~<r»> 




Figure 7: The distribution of gluons ordered in rapidity, as in our basic LO MC. 



The reason why POWHEG technique is complicated in case of the angular ordering is illus- 
trated in Fig. 7. We show there the distribution of gluons ordered in rapidity, starting from the 
gluon with the maximum rapidity, the closest to hard process. The gluon distribution with the 
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highest rapidity ^ ~ i^max (•/=!) has a ridge extending towards the soft region. Notice that, when 
the IR cut-off e — )• in (1 — z) < £, the width of this ridge also goes to zero. Consequently, the 
gluon with the highest ^ is unable to reproduce the gluon distribution in the NLO corner, close to 
hard process. This is why in this case POWHEG requires truncated and vetoed showers, which are 
not needed in our method. 

5. Summary and outlook 

A new method of adding the QCD NLO corrections to the hard process in the initial state 
Monte Carlo parton shower is tested numerically showing that the basic concept of the new method- 
ology works correctly in the numerical environment of a Monte Carlo parton shower. The differ- 
ences with the well established methods of MC@NLO and POWHEG are briefly discussed. Also, 
variants of the new method with better efficiency in terms of CPU time are proposed. 
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